Search for: All records

We present an implementation of the relativistic ionization-potential (IP) equation-of-motion coupled-cluster (EOMCC) with up to 3-hole–2-particle (3h2p) excitations that makes use of the molecular mean-field exact two-component framework and the full Dirac–Coulomb–Breit Hamiltonian. The closed-shell nature of the reference state in an X2C-IP-EOMCC calculation allows for accurate predictions of spin–orbit splittings in open-shell molecules without breaking degeneracies, as would occur in an excitation-energy EOMCC calculation carried out directly on an unrestricted open-shell reference. We apply X2C-IP-EOMCC to the ground and first excited states of the HCCX+ (X = Cl, Br, I) cations, where it is demonstrated that a large basis set (i.e., quadruple-zeta quality) and 3h2p correlation effects are necessary for accurate absolute energetics. The maximum error in calculated adiabatic IPs is on the order of 0.1 eV, whereas spin–orbit splittings themselves are accurate to ≈0.01 eV, as compared to experimentally obtained values. 

Recent theoretical studies have explored how ultra-strong light–matter coupling can be used as a handle to control chemical transformations. Ab initio cavity quantum electrodynamics calculations demonstrate that large changes to reaction energies or barrier heights can be realized by coupling electronic degrees of freedom to vacuum fluctuations associated with an optical cavity mode, provided that large enough coupling strengths can be achieved. In many cases, the cavity effects display a pronounced orientational dependence. Here, we highlight the critical role that geometry relaxation can play in such studies. As an example, we consider a recent work [Pavošević et al., Nat. Commun. 14, 2766 (2023)] that explored the influence of an optical cavity on Diels–Alder cycloaddition reactions and reported large changes to reaction enthalpies and barrier heights, as well as the observation that changes in orientation can inhibit the reaction or select for one reaction product or another. Those calculations used fixed molecular geometries optimized in the absence of the cavity and fixed relative orientations of the molecules and the cavity mode polarization axis. Here, we show that when given a chance to relax in the presence of the cavity, the molecular species reorient in a way that eliminates the orientational dependence. Moreover, in this case, we find that qualitatively different conclusions regarding the impact of the cavity on the thermodynamics of the reaction can be drawn from calculations that consider relaxed vs unrelaxed molecular structures. 

Polariton chemistry exploits the strong interaction between quantized excitations in molecules and quantized photon states in optical cavities to affect chemical reactivity. Molecular polaritons have been experimentally realized by the coupling of electronic, vibrational, and rovibrational transitions to photon modes, which has spurred a tremendous theoretical effort to model and explain how polariton formation can influence chemistry. This tutorial review focuses on computational approaches for the electronic strong coupling problem through the combination of familiar techniques from ab initio electronic structure theory and cavity quantum electrodynamics, toward the goal of supplying predictive theories for polariton chemistry. Our aim is to emphasize the relevant theoretical details with enough clarity for newcomers to the field to follow, and to present simple and practical code examples to catalyze further development work.